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Abstract 

The direct observation of gravitational waves will provide a unique tool for probing the dynamical 
properties of highly compact astrophysical objects, mapping ultra-relativistic regions of space-time, 
and testing Einstein's general theory of relativity. LISA (Laser Interferometer Space Antenna), 
a joint NASA-ESA mission to be launched in the next decade, will perform these scientific tasks 
by detecting and studying low-frequency cosmic gravitational waves through their influence on 
the phases of six modulated laser beams exchanged between three remote spacecraft. By directly 
measuring the polarization components of the waves LISA will detect, we will be able to test 
Einstein's theory of relativity with good sensitivity. Since a gravitational wave signal predicted by 
the most general relativistic metric theory of gravity accounts for six polarization modes (the usual 
two Einstein's tensor polarizations as well as two vector and two scalar wave components), we have 
derived the LISA Time-Delay Interferometric responses and estimated their sensitivities to vector- 
and scalar-type waves. We find that (i) at frequencies larger then roughly the inverse of the one- 
way light time (~ 6 x 10^^ Hz.) LISA is more than ten times sensitive to scalar-longitudinal and 
vector signals than to tensor and scalar-transverse waves, and (ii) in the low part of its frequency 
band is equally sensitive to tensor and vector waves and somewhat less sensitive to scalar signals. 
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I. INTRODUCTION 



The direct detection of gravitational waves will represent one of the greatest triumphs 
of experimental physics of this century, and provide us with a new observational tool for 
obtaining better and deeper understanding about their sources. 

Several experimental efforts have been underway for years, both on the ground and in 
space [IHS], and only recently kilometer-size ground-based interferometers have been able 
to identify the most stringent upper-limits to date for the amplitudes of the radiation ex- 
pected from several classes of sources. Although an unambiguous detection has not been 
declared with present-generation instruments, next-generation Earth-based interferometers 
and pulsar-timing experiments, as well as the LISA (Laser Interferometric Space Antenna) 
mission are expected to achieve this goal. 

LISA, jointly proposed to the National Aeronautics and Space Administration (NASA) 
and the European Space Agency (ESA), is expected to be flown sometimes in the next 
decade. Its goal is to detect and study gravitational waves (GW) in the millihertz frequency 
band. It will use coherent laser beams exchanged between three identical spacecraft forming 
a giant (almost) equilateral triangle of side 5 x 10^ kilometers. By monitoring the relative 
frequency changes of the light beams exchanged between the spacecraft, it will extract the 
information about the gravitational waves it will observe at unprecedented sensitivities [?]. 
The astrophysical sources that LISA is expected to observe within its operational frequency 
band (10""^ — 1 Hz) will be very large in number, including galactic and extra-galactic 
coalescing binary systems containing white dwarfs and neutron stars, extra-galactic super- 
massive black-hole coalescing binaries, and a stochastic gravitational wave background from 
the early universe [H [9] . 

The first unambiguous detection of a gravitational wave signal, whether performed on 
the ground or in space, will also allow us to test Einstein's general theory of relativity (GR) 
by measuring the polarization components of the detected signals [101 [H]. Among all the 
proposed relativistic metric theories of gravity [12], GR is the most restrictive, allowing for 
only two of possible six different polarizations [13] . By asserting that the spin-2 ( "tensor" ) 
polarizations are the only polarization components observed, we would make a powerful 
proof of the validity of Einstein's theory of relativity, while a clear observation of some other 
polarization modes would disqualify it. Corroboration of polarization measurements with 
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estimates of the propagation speed of the observed gravitational wave signal will provide 
further insight into the nature of the observed radiation and result into the determination 
of the mass of the graviton. It should be emphasized, however, that a propagation-speed 
measurement alone consistent with a value equal to the speed of light, would not automati- 
cally rule-out other relativistic metric theories of gravity. This is because waves with helicity 
s = can also propagate at light-speed [H]. Gravitational waves with scalar polarization 
are predicted by the most common generalizations of GR such as scalar-tensor theories. 
Besides the classical example of Brans-Dicke theory [T3] , scalar-tensor theories result in the 
low-energy limit of string theory (see, e.g., [IS]). Modifications of the Einstein-Hilbert ac- 
tion, which consider generic functions of the Ricci scalar in the Lagrangian (f(R) theories), 
also predict "scalar" gravitational waves [17] . Vector modes, on the other hand, can appear 
in the so-called "quadratic gravity" formulations [17] , and in the context of theories in which 
the graviton has a finite mass such as the Visser theory |18j . 

LISA will not be able to distinguish the propagation speeds of scalar (helicity s = 0) 
and vector (helicity s = ±1) polarizations from the speed of light (a result following, as we 
shall see below, from a combination of the existing stringent upper limits on the mass of 
the graviton [T9H2T] and the LISA observational bandwidth). However, it should be able 
to assess the polarizations of the observed gravitational wave signals. Since the accuracy 
by which LISA will distinguish one polarization from another will depend on their signal- 
to-noise ratios [22], in this paper we estimate the LISA sensitivities to scalar- and vector- 
polarized wave. 

The paper is organized as follows. In Section (|ll]) we derive the one-way Doppler response 
to a gravitational wave signal characterized by six polarizations (2 "tensor" (helicity s = ±2), 
2 "vector" (helicity s = ±1), and 2 "scalar" (helicity s = 0)). Since the resulting expression 
is equal in form to the one previously derived by Estabrook and Wahlquist [23] for tensor 
waves (i.e. for waves predicted by General Relativity), we conclude that the responses of the 
various time-delay interferometric (TDI) combinations are also identical in form to those 
previously derived within the framework of GR [21] . Although the derivation of the response 
function of a Michelson interferometer to non-tensor polarization modes has been considered 
in previous publications [TOj |25] , the expression presented there was correct in the so called 
"long- wavelength-limit" , i.e. when the wavelength of the GW is much larger than the size of 
the detector. Since LISA will be sensitive, over most of its observational frequency band, to 
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GWs of wavelength shorter than its hnear size, we have derived the expressions of the LISA 
TDI responses to vector- and scalar- waves that are valid for any wavelength. In Section 



(III) we then compute the LISA sensitivities to vector and scalar waves. We find that (i) 
at frequencies larger then the inverse of the one-way light time 6 x 10^^ Hz) LISA is 
ten-times more sensitive to scalar-longitudinal and vector signals than to tensor and scalar- 
transverse waves, while (ii) it is equally sensitive to tensor and vector waves, and somewhat 
less sensitive to scalar signals in the low part of its frequency band. Although both these 
results might seem surprising at first, we show that they are consequence of the physical and 
geometrical properties of the vector and scalar waves and how they affect the frequency of 



the light beams exchanged by the three LISA spacecraft. Finally in Section (IV) we present 
a summary of the paper and our concluding remarks. Throughout the paper we will be 
using natural units {c = G = h = l) except where mentioned otherwise. 



II. DERIVATION OF THE ONE-WAY DOPPLER RESPONSE 

In what follows we present the derivation of the "one-way" Doppler response to a gravi- 
tational wave signal predicted by the most general relativistic metric theory of gravity [T3] . 
Although the result has already appeared in the literature PH], our derivation results into an 
expression that is compact and identical in form to that first obtained by Wahlquist in 
the case of plane gravitational waves predicted by GR. This of course simplifies significantly 
the derivation of the LISA TDI responses as they turn out to be equal in form to those 
given in pi] . As in [231 [2Z] , our derivation is general and does not rely on any assumptions 
about the size of the wavelength of the radiation relative to the size of the detector. 

Let us consider a space-time with the following metric 

ds^ = —dt^ + {5ij + hijiyt — z))dx'^dx\ \hij\ « 1, (1) 

where the usual sum convention over repeated indices is assumed, Latin indices go from 
1 to 3, and a plane- gravitational wave has been assumed, without loss of generality, to 
propagate along an arbitrary +z direction. In Eq. ([T| we have allowed the wave to travel 
at a finite speed (group velocity) f < 1 to account for a possible non-zero mass of the 
graviton. Working in the context of GR, it is well know that hij has two degrees of freedom 
representing gravitational waves (GWs) with helicity s = ±2, and v = 1. On the other 
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hand, alternative relativistic metric theories of gravity allow for GWs with up to six degrees 
of freedom [13]. Therefore, in the most general case, hij can be represented in terms of six 
components (with corresponding six metric amplitudes) in the following form [13] 

6 

hij{vt- z) = Y^ h^r) {vt~z) , (2) 



r=l 



(t) 

where e\j are the six polarization tensors associated with the six waveforms of the grav- 
itational wave signal. If we introduce a set of Cartesian orthogonal coordinates (x, y, z) 
associated with the wave, in which (x, y) are in the plane of the wave and z is along the 
direction of propagation of the wave and orthogonal to the (x,?/) plane, the above six po- 
larization tensors assume the following matricial form 
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From the above expressions it is easy to verify that the tensors e*^") , a = 1, 6 are linearly 
independent and form an orthogonal basis. In our notation we have labeled 4 and 5 the 
usual + and x polarizations respectively (the tensor polarization waveforms). The vector 
polarizations (s = ±1) were labeled as 2 and 3, and finally the two scalar modes (s = 0) 
have been denoted with the labels 1 (for the longitudinal scalar polarization) and 6 (for the 
transversal scalar mode). 

Following [231 [28] we may notice that the space-time described by the metric of Eq. ([T]) 
allows for the following three Killing vectors 

<5f 



^(1) 



^(2) 



^(3) 



5f + 



(4) 



Generally speaking, in the weak field regime of a generic metric theory of gravity, GWs can 
travel with f < 1. In that case the linearized field equations are Klein-Gordon-type, and the 
resulting group velocity is determined by the mass of the graviton, m. Remembering that 
V = dtu/dk, and that the dispersion relation is equal to k = y/co"^ — m^, we find the following 
expression for v in terms of m and u 

- ■ («) 

Since the operational frequency band of LISA will be within the range (10"'* — 1 Hz), by 
assuming presently known upper-limits for the mass of the graviton from Eq. ^ we obtain 
the resulting values for the group velocity of these waves. After restoring physical units and 
taking m < 10^^^ g (the most stringent constraint to date obtained by requiring the derived 
dynamical properties of a galactic disk to be consistent with observations [21]), at 10~^ Hz 
we find a value for the group velocity v whose fractional difference, A^,, from the speed of 
light is equal to: = |i; — c|/c ~ 10~^^. A less stringent value for the mass of the graviton 
equal to m < 7.68 x 10^^^ g (obtained from solar-system dynamics observations [E]), results 
into a At, ^ 10^®. These considerations imply that, no matter whether we are conservative 
or not in our assumption about a likely upper limit for the graviton mass, LISA will not 
be able to resolve the propagation speeds of the different polarization components of the 
detected GWs by resolving the time-separations of their imprints in the TDI combinations 
[21]. For this reason, from now on, we will assume f = 1 in natural units, and rewrite the 
third Killing vector 1^(3) as K^^^^ = 6^ + 6^. 

Let us now consider the unit vector n along the Doppler link and oriented from spacecraft 
1 to 2 (see Fig. [T]). If we denote with (6',0) the usual polar angles, h assumes the following 
familiar form 

n'' = sin 9 cos 4> + sin 9 sin (p Sy + cos 9 51 . (6) 

Since the most general Killing vector of our metric can be written as a linear combination 
of the three Killing vectors above 

RP = aiSP + a2SP + as{SP + (5f) , (7) 

with ai, 02, fla constants, by comparing Eq. ^ and Eq. ([t]) we note that by taking 
(01,02,03) = C n, Eq. ^ becomes 

= C[nP + k-h6P] , (8) 
6 



FIG. 1: A laser beam of nominal frequency i'q is transmitted from spacecraft 1 to spacecraft 2 and 
simultaneously another beam from spacecraft 2 is transmitted back to 1. The gravitational wave 
train propagates along the z direction, and the two polar angles {9, (p) describe the direction of 
propagation of the laser beams relative to the wave. See text for a complete description. 

where C is an arbitrary constant, and k is the unit vector along the direction of propagation 
of the wave (see Figure [T]). 

If we now consider the 4-momentum of a photon transmitted by spacecraft 1, its analytic 
expression can be written in terms of the metric perturbation at spacecraft 1 in the following 
way [231 EH] 

Pp = i^o(^-Sl + np + hp^n^^ , (9) 

where it is easy to see that the condition PpP^ = is fulfilled to first order in the metric 
perturbation hij. Since PpK^ = constant along the photon world line [2H], we obtain 
the following relationship between the frequency of the photon emitted at spacecraft 1, z/q, 
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and that received at spacecraft 2, v' 



P,KP = P'pK'P, 



(10) 



1^0 (^-Sl + np+ ^ j {n^ + k- h6?) = v' + n; + \h:p^n^ [n" + k ■ n'bl) . (11) 

If we now rewrite n!'' = nP + 6n'', we may notice that, to first order in hij, n^dnp = because 



' nP 



n„n' 



npU^ = 1. This resuh, together with Eq.(ll) above, allows us to derive the following 



expression for the ratio of the two frequencies u' and z/q 



u' 1 — k ■ h + ^n^hp^n^ 
1 — k ■ n + ^nPh'p^n^ 



(12) 



Finally, expanding the right-hand-side of the last equation to first order in hp^ we get the 
resulting expression for the one-way Doppler response, y 
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where y(t) = {^'if) — vq)Ivq, and \l/(t) is equal to 



(13) 
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2[1 - {k-hy\ 

By explicitly showing the time dependence of the various terms, ?/()!:) can be rewritten in the 
following form [2] 



'\^k-n) 



^{t-L)--^{t-k-hV) 



(15) 



where L is the separation between the two spacecraft. Note that, in order to obtain the 
above expression, we have only assumed the time-components of the metric perturbation, 
hpt to be equal to zero [13]. 

The expression for the one-way Doppler response measured on board spacecraft 1 at time 



t can be obtained from Eq. ( 15 ) by changing fi — )■ — n, and further delaying the waveforms 



hj k ■ n L. The resulting one-way Doppler response, y'it), is equal to 

y'{t) = {l-k-n) h{t -{l + k- h)L) - ^{t) 



(16) 



Since the above expressions of the one-way Doppler responses are identical in form to 
those valid for tensor waves [27], it follows that the various TDI combinations of the LISA 
six inter-spacecraft one-way Doppler measurements will also be identical in form to those 
derived in [21]. For these reasons they will not be given here, and we refer the reader to 
El [291 130] for more details. 
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III. LISA SENSITIVITIES 



In this section we will compute the LISA TDI sensitivities to vector and scalar waves. 
Although the LISA sensitivities to tensor waves have already been presented in an earlier 
publication [31j, for sake of comparison we will include them in the sensitivity plots presented 
in this section. We will specialize our calculations to the equilateral-triangle configuration: 
armlength 1 = armlength 2 = armlength 3 = L (L ~ 16.7 light seconds) since the LISA 
arm lengths will differ by at most a few percent, and any corrections to our results will be 
to this level of accuracy |211 EHl EI] • 

The LISA sensitivity to tensor GWs has been traditionally taken to be equal to (on 
average over the sky and polarization states) the strength of a sinusoidal GW required to 
achieve a signal-to-noise ratio of 5 in a one-year integration time, as a function of Fourier 
frequency. Although in the case of vector waves the average over the polarization states can 
be performed by implementing the same procedure described in [2H [31] for tensor waves, 
in general this can not be done for scalar waves. This is because the two scalar fields are 
mutually orthogonal (and independent), as one is purely longitudinal and the other purely 
transverse to the direction of propagation of the wave. 

In order to compute the LISA sensitivities we will use the following expressions for the 
power spectra of the noises affecting the X, a, (, E, P, and U combinations [31] (see Figure 
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FIG. 2: Noise spectra in the X, a, C) ^'^^ U time-delay interferometric combinations 

accounting for the nominal proof- mass and optical path noises. The varying depths of the minima 
in the high-frequency ranges of X, E, and P is an artifact of numerically calculating these functions 
at discrete frequencies 

where S^^{f) = 2.5 x 10~^^ [f/lHz]~'^ Hz~^ is the spectrum of the relative frequency 
fluctuations due to each proof mass, and S°P{f) = 1.8 x 10-^'^[f /IHzf Hz"^ is the spectrum 
of optical path (mainly shot and beam pointing) noise. Both these noises can be regarded 
as the main limiting noise sources for LISA [3 [31] . 

Gravitational wave sensitivity is the wave amplitude required to achieve a given signal-to- 
noise ratio. We calculate it in the conventional way, requiring a signal-to-noise ratio of 5 in a 
one year integration time: 5 ^/ Sk{f) B/ (root-mean-squared gravitational wave response for 
data combination k), where k is X, a, (, E, P, U, and Sk is the total noise power spectrum 
for TDI combination k. The bandwidth, B, was taken to be equal to one cycle/year (i.e. 
3.17 X 10-^ Hz). 

We have assumed the vector-waves to be elliptically polarized and monochromatic, with 
their wave functions, {h^'^\ h^^''), written in terms of a nominal wave amplitude, H, and the 
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two Poincare parameters, F), in the following way 



/;,(2)(t) = H sin(r) sin(u;t + <l>) , 
h^^\t) = H cos(r) sin(wt) . 



(23) 
(24) 



For scalar signals instead, the two wave functions, {h^^\ h^^^), have been treated as indepen- 
dent and we calculated the TDI sensitivities to each of these two polarizations. For both 
vector and scalar signals we averaged over source direction by assuming uniform distribution 
of the sources over the celestial sphere; in the case of vector signals we also averaged over 
elliptical polarization states uniformly distributed on the Poincare sphere for each source 
direction. The averaging was done via Monte Carlo integration with 4000 source posi- 
tion/polarization state pairs per Fourier frequency bin and 7000 Fourier bins across the 



Figure ^ shows the root-mean-squared (r.m.s.) responses of the TDI combinations (a) 
X (unequal-arm Michelson), (b) a (Sagnac), (c) ( (symmetrized Sagnac), (d) E (monitor), 
(e) P (beacon), and (f) U (relay) to tensor (already derived in [24j and shown here for 
comparison), vector and scalar gravitational waves. In the high-part of the LISA frequency 
band we may notice that the r.m.s. responses to vector and scalar-longitudinal waves are 
significantly larger than those to tensor and scalar-transverse signals. In particular, the 
scalar-longitudinal r.m.s. response grows with the frequency at a much faster rate than the 
others. 

In order to physically understand this effect, let us compare a one-way Doppler response 
(that measured on board spacecraft 1, for instance) to a "pulse" tensor wave against that due 
to a scalar- longitudinal pulse wave. A tensor signal propagating orthogonally to the light 
beam (direction for which the one-way Doppler response can reach its maximum magnitude 
in this case) will only interact with the light for the brief instance its wavefront crosses the 
light beam. On the other hand, if a scalar-longitudinal wave propagates along the direction 
between the two spacecraft (over which the Doppler response will achieve its maximum in 
this case) the frequency of the light will be affected by the gravitational wave for the entire 
time L it takes the wave to propagate from one spacecraft to the other, resulting into an 
amplification of the frequency shift when the wavelength of the wave is shorter than the 
inter-spacecraft distance, L. The above considerations become apparent by considering in 
both cases the one-way Doppler response y' in the Fourier domain. For a tensor signal, the 
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modulus-squared of the Fourier transform of Eq. (16) with k ■ n ^ becomes equal to 



2 



sm\7rfL) , 



(25) 



fc-n— ^0 



with h = /i(^) cos(20) + /i(^) sin(20). In the case of a scalar-longitudinal signal instead, it is 
easy to show that, in the limit of A; ■ n — ?■ — 1 the modulus-squared of the Fourier transform 



of Eq. (16) becomes equal to 



y'if)k 



k-n- 
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(26) 



From the above two expressions we may conclude that, at frequencies larger than the inverse 
of the one-way-light-time and for tensor and scalar-longitudinal waves of comparable ampli- 
tudes, the maximum one-way Doppler response to a scalar-longitudinal signal will be larger 
than the corresponding maximum response to a tensor wave by roughly a factor of vr/L. 
This example implies that the r.m.s. of the TDI responses to scalar-longitudinal signals will 
be larger than those to tensor waves in the high-frequency region of the LISA band. Similar 
considerations can be made for understanding the differences between the r.m.s. responses 
to tensor and vector waves. 

In the low-frequency limit (for / < 5 x 10^'^ Hz) the tensor and vector r.m.s. responses 
coincide, while those for the two scalars also coincide with each other but are smaller by 
a factor of about 2 than those for tensor and vector waves. This is because the two scalar 
waves, being mutually orthogonal, have been treated as independent rather than elliptically 
polarized like the vector and tensor waves. 

In Figure ^ we then plot the corresponding sensitivities of the TDI combinations to the 
tensor, vector and the two scalar polarization components. The characteristic behavior of 
the r.m.s. responses to scalar-longitudinal and vector signals shown in figure ^ folds into 
the plots presented here. At high-frequencies the sensitivity to scalar-longitudinal waves is 
significantly better than that to tensor, vector, and scalar-transverse waves. At 1 Hz, for 
instance, the sensitivity of scalar-longitudinal signals is about a factor of fL ^ 17 better 
than that to tensor waves, while it is only a factor of 3 better than the sensitivity to vector 
signals. 

Another interesting feature noticeable in Figure ^ is the lack of sensitivity of the Sagnac 
combinations (b) a and (c) ( to scalar-transverse waves at frequencies equal to integer- 
multiples of the inverse of the one- way- light-time. We have verified this result analytically 
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and found that indeed, at these frequencies, the Sagnac responses are identically equal to 
zero independently of the direction of propagation of the signal. 



IV. SUMMARY AND CONCLUSIONS 

The main results of our work have been that (i) LISA is more sensitive to scalar- 
longitudinal and vector signals than to tensor waves in the high-part of its frequency band, 
and (ii) at low-frequencies its sensitivities to tensor and vector signals are equal and some- 
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what better than those to scalar waves. We have also found that the LISA TDI Doppler 
responses to scalar-longitudinal waves propagating along any of the three LISA arms will 
experience an amplification proportional to the arm-length. These results, together with 
the LISA capability for constructing three independent TDI combinations in the high part 
of its observational frequency band, should provide LISA with the capability for assess- 
ing the polarization of the waves it will detect. This will be the topic of our forthcoming 
investigation. 
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FIG. 4: Sensitivities of the (a) X, (b) a, (c) (d) E, (e) P, and (f) U TDI combinations to 
gravitational waves with tensor (s = ±2), vector (s = ±1), and scalar (s = 0) components. 
Consistently with Figure ([s]) we may notice how more sensitive the various TDI combinations 
are to vector and scalar-longitudinal signals than to tensor and scalar-transverse waves in the 
high-frequency part of the LISA frequency band. See text for more details. 
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